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Abstract 

Recently several mathematical models of a solar chimney power plant 
were derived, studied for a variety of boundary conditions, and compared 
against CFD calculations. The importance of these analyses is about the 
accuracy of the derived pressure drop and output power equation for solar 
chimney power plant systems (SCPPS). We examine the assumptions under¬ 
lying the derivation and present reasons to believe that some of the derived 
equations, specihcally the power equation in this model, may require a cor¬ 
rection to be applicable in more realistic conditions. The analytical resutls 
are compared against the available experimental data from the Manzanares 
power plant. 
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Nomenclature 

Variables 

A cross-sectional area 

Aj. cross-sectional area of the collector ground 

g acceleration due to gravity 

h height 

fh air mass flow rate 
p pressure 

W flow power 
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q heat transfer per unit mass 

q" heat flux 

R gas constant 

T temperature 

p density 

u velocity 

Cp specific heat capacity 

Subscripts 
i inlet 

o outlet 

c collector 

t tower 

turb turbine 

atm atmospheric 

Abbreviations 
LHS left hand side 
RHS right hand side 


1. Introduction 

Although the idea of the solar chimney power plant (SCPP) can be traced 
to the early 20th century, practical investigations on solar power plant sys¬ 
tems started in the late 1970s, around the time of conception and construction 
of the hrst prototype in Manzanares, Spain. This solar power plant operated 
between 1982 and 1989 and the generated electric power was used in the local 
electric network [1, 2], 

The basic SCPP concept (Fig. 1) demonstrated in that facility is fairly 
striaghtforward. Sunshine heats the air beneath a transparent roofed col¬ 
lector structure surrounding the central base of a tall chimney tower. The 
hot air produces an updraft flow in the chimney. The energy of this up¬ 
draft flow is harvested with a turbine in the chimney, producing electricity. 
Experiments with the prototype proved the concept to be viable, and pro¬ 
vided data used by a variety of later researchers. A major motivation for 
subsequent studies lay in the need for reliable modeling of the operation of a 
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large-scale power plant. The Manzanares prototype had a 200 m tall chim¬ 
ney and a 40,000 collector area. Proposals for economically competitive 
SCPP facilities usually feature chimneys on the scale of 1 km and collectors 
with multiple square kilometer areas. 



Figure 1: Schematic of SCPP with the applied variables and subscripts in the present 
analysis 

Padki and Sherif [3] used the results from the Manzanares prototype to 
extrapolate the data to large scale models for SCPP. In 1991, Yan et al. [4] 
developed an SCPP model using a practical correlation. They introduced 
equations including air velocity, air flow rate, output power, and thermofluid 
efficiency. Von Backstrom and Fluri conducted a numerical study to de¬ 
termine the optimum ratio of pressure drop of the turbine as a fraction of 
the available pressure difference required to achieve the maximum power [5] . 
They noted that this ratio might lead to overestimating the flow passage 
in the plant and also designing a turbine without a sufficient stall margin. 
In other recent works, the SCPP concept involving an inflatable tower was 
examined, with all parts of the power plant modeled numerically [6, 7, 8]. A 
small-scale inflatable tower was fabricated for validation of these results, and 
code calibration was performed using the newly available experimental data 
[9, 10]. 

To hnd the maximum power, different atmospheric pressure and temper¬ 
ature boundary conditions were applied for various tower heights and atmo¬ 
spheric lapse rates [11]. Theoretical analysis to study the effect of pressure 


3 













drop in the SCPP turbine was performed by Koonsrisuk et al. [12], The 
optimal pressure drop ratio was found numerically and analytically by Gue 
et ah, around 0.9 for the Manzanares prototype. This inverstigation can be 
applied as an initial estimation for various SCPP turbines [13]. Tayebi et al. 
modeled and simulated the SCPP with a curved conjunction between tower 
and collector for different Rayleigh numbers [14]. 

Earlier modeling efforts [7] showed a keen sensitivity of the predictions 
of SCPP output to boundary conditions, in particular, pressure. Numerical 
simulations require careful validation and verihcation, and for that, analytical 
models are indispensable. A theoretical model was recently developed [15] to 
model the combined performance of the solar collector, chimney, and turbine. 
Here we will examine some of the assumptions and derivations in this model 
and present an alternative formulation for the energy equation. 

2. Derivation of equations 

2.1. Collector 

To derive the equations, we start from the collector. It is assumed that 
the flow through the collector is one-dimensional, steady-state, and com¬ 
pressible. Let us disregard the friction and assume the total heat from the 
solar irradiation is absorbed within the air hlling the collector. For this one 
dimensional axisymmetric compressible flow analysis, the mass conservation 
satishes: 


dA dp du 

——I-1- 

A p u 


0 (Continuity) 


( 1 ) 


Here A is the cross-sectional area of the collector that air goes through ~ 
A = 2'Krhc and dA = 2'Krdhc. 

Momentum equation is as follows [16]: 

The momentum equation can be obtained from: 



dr dr 


( 2 ) 


dp + pudu = 0 (Momentum) 
The energy equation and the equation of states are 

pude pdu kd‘^T 

+ ^ + = 0 


dr dr dx'^ 


( 3 ) 


( 4 ) 
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P = p{p,T),e = e{p,T) (5) 

The last two equations, (5), represent the thermal and caloric equations 
of state [17]. Consider the energy balance equation and the equation of state 
as follows: 


CpdT — dq + udu = 0 (Energy) 


( 6 ) 


+ f )-0 (State) 
P \P T ' 


(7) 


To hnd dp we can apply Eq. (3) and substitute du/u from the continuity 
equation, Eq. (1). 


dp du dp dA 
P 


pu‘^ 


dp = pu [— + ^ 


u 


2 I dp ^ dA 


( 8 ) 


(9) 


From the equation of state we can hnd dp/p and substitute in Eq. (9), 

( 10 ) 


dp dp dT 
p p T 


2 / dp dT ^ dA 


dp = pn + ^ 


( 11 ) 


We can rewrite Eq. (11) as a function of T, A, u^p, p , where m = pAu\ 

( 12 ) 


^ p‘^v?A^ ^ dp dT ^ dA^ 


p \A‘^p A^T 713 J 

Also by substitution dT from the energy equation on the base of dq, Cp 
and u, we obtain 


m 


dp = — ypT 

p 


f dA dq — udu dp 


+ 


(13) 


VA3 A'^Tcp ' A^p^ 

For consistency with previous analyses, let us rewrite dq on the basis of 
heat hux per mass how rate— dq = q"dA^/rn where q has the units of J/kg. 
Here Ar = vrr^, therefore dAr = 27rrdr. Note that A = 27rrhr, where hr is 
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the collector height (roof height) that was assumed to be proportional to r - 
he = ar, where a is a constant. By substituting Ar, dq and A in the second 
term on the RHS, we obtain 


• 2 

m 


dp = — {^ 
P 


q"{2'Kr)dr 


+ 


udu 

A^CpT 


+ 


dp 


(14) 


v43 m{2'Kr‘^afTcp A^CpT A^p ^ 

We can rewrite equation (14) and substitute udu of the third term on the 


RHS by applying momentum equation (3), udu = —dp/p. 


• 2 

m 


dp='—(^- 
P 


q''dr 


dp 


+ 


dp 


(15) 


A^ 2'Kmr^a?CpT A^pCpT A^p ^ 

Then we can substitute p from the equation of state, p = pRT and rewrite 
the above equation to hnd dp on the LHS. 


• 2 

m 

'dA 

1- 

1 

m 


01 

p 


2'Kmr^a^CpT 

A^p^T 

KR 

CpJ. 


Note that 


m = pHu, 


(16) 


So we can rewrite (16) as 


' 2 

m 

'dA 



[i 


01 

p 


2TTrnr^o?‘CpT 


T 

u 

Cp/. 


(17) 


Equations (16) and (17) are the exact solutions for dp for the one-dimensional 
frictionless analysis of the collector. Since our fluid is air we can estimate Cp 
and rewrite Eq. (17). 


TV? f dA q"dr \ / 2.494n^\ ^ 

On the basis of mass flow rate, Eq. (18) can be written as: 


(18) 


/dA q''dr \ / 2.494m^\ ^ 

^ “ y “ 27rmr^a^CpT ) ~ Tp^A^ ) 

The third term of the RHS of Eq. (19) was ignored [15] which can be 
correct when density is constant. Cp, q" and T are considered approximately 
constant as well. Therefore by integrating between the inlet and outlet of 
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the collector without the last term of the RHS, pressure difference can be 
derived. 



mq''dr \ 
2T[r^a?pCpT ) 


r • 9 

m 1 


1 qm 


UL AlJ 

Ana^pCpT 




( 20 ) 

( 21 ) 


2.2. Tower 

The air flow in the tower(chimney) is considered as an adiabatic friction¬ 
less flow. The conservation equations for the one-dimensional steady state 
flow in variable-area tower are as follows: 


dp du dA r. . . N 

-1-1—— = 0 (Continuity) 

p u A 

dP 

- h gdz -|- udu = 0 (Momentum) 

CpdT + udu -|- gdz = 0 (Energy) 

dp_dp_^^^ (State) 

p p T 

By following the same trend to hnd dp we get 


dp = 


—pgdz + 


m^dA 

pA^ 


dp dT 
+ P“ — - ^ 


P 


T 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


By applying the energy equation and substitution dT = {—gdz — udu)/cp, 
we can rewrite the above equation as 


dp = 


—pgdz + 


rh^dA 

pA^ 


+ pu -h 


dp gdz -|- udu 


p 


CpT 


Here dp = —p{udu + gdz), then we get 
dp = 


, m^dA 

-pgdz + —— + puU — 


pA^ 

Above equation can be solved for dp, 


2 fdp _ dp \ 

p pCpT) 


(27) 


(28) 
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-1 


dp = 


-pgdz + 


rfi^dA 

pA^ 


u 


T \R 


1-^ 


(29) 


Also by considering the material properties of air the same way we did 
for the collector part, 


dp ~ 


—pgdz + 


rn?dA 

pA^ 



2.494n2]-^ 

T 


(30) 


The above eqnation is the exact closed form solntion of dp at any point on 
the base of variable p, T and q". We ignore the last term on the RHS by 
having a constant density [15] and integrate between the inlet and ontlet 
tower area to hnd the pressnre difference of the chimney as, 



(^pgdz + 


m^dA\ 

~w) 


(31) 


Pt,i - Pt,o + pght + 
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m 




(32) 


To calcnlate the ontpnt power, we can dehne the power on the basis of 
the pressnre difference at the tnrbine ~ where it is normally ntilized at the 
ontlet of the collector and inlet of the tower. 


W ~ 


m{pc,o - Pt,i) 

Pturb 


(33) 


Let Pturb = {Pc,o + Pt,i)/‘^ and snbstitnte eqnations Pc,o and pt^t from (21) 
and (32). Hence for the flow power by assnming = Pt,o + pght, we have 


W 


m 


(,Pc,o T Pt,i)/‘^ 


-m 


2p 


^lo 


1 \ q''m ( 1 
^c,i/ ^ 47ra2pCpT Vr2^ 



wA ( 1 



( 34 ) 


For area the following eqnations are nsed, where b and c are arbitrary 
positive real constants. 


^c,i ~ ^t,o ~ ( 35 ) 

The simplihed form of eqnation (34) by applying the area correlations is. 
























w ~ 


i,Pc,o Pt,i)/‘^ 2p 


m 


• 2 
—m 



(36) 


Assume p is constant before and after turbine, therefore 


Pturb Pc,o Pt,i P 


Koonsrisuk et al. derived an equation in which the second term was 
neglected in comparison with the hrst term on the RHS of Eq. (36). However, 
Eq. (37) shows the derived power equation by them at the end is likely to 
exceed the expected amount by a factor of two. 



(37) 


To evaluate the derived anlaytical solution for the output power of SCCP, 
the available experimental data from Manzanares prototype was applied and 
extracted. The measured updraft velocity for 25 hours Manzanares power 
plant operation is imposed to the analytical solution and the analytical power 
compared against the experimental outpower from the turbine (Fig. 2). 



Figure 2: Analytical power results against measurement from Manzanares: updraft veloc¬ 
ity and power output for a typical day 
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3. Conclusion 


We presented considerations supporting our contention that a known an¬ 
alytical model predicting a solar chimney power plant performance may have 
a power equation which is off by a factor of two. Careful derivation of the 
models is very important, especially for the specihc area of interest related to 
solar-chimney power plants, where numerical model scalability is a key issue, 
and few experimental results are available for validation. During the verifi¬ 
cation and validation process, the modeler must ask two questions: ”Am I 
modeling the physics correctly?” and ”Am I modeling the correct physics?” 
Comparison with analytical models is important for answering both of these 
questions, and the only way to have them well-posed is to have correct physics 
in the analytics. 

References 

References 

[1] Haaf, W., K. Friedrich, G. Mayr, and J. Schlaich. ’’Solar chimneys part 
I: principle and construction of the pilot plant in Manzanares.” Interna¬ 
tional Journal of Solar Energy 2, no. 1 (1983): 3-20. 

[2] Haaf, W. ’’Solar chimneys: part ii: preliminary test results from the 
Manzanares pilot plant.” International Journal of Sustainable Energy 2, 
no. 2 (1984): 141-161. 

[3] Padki, M. M., and S. A. Sherif. ’’Solar chimney for medium-to-large scale 
power generation.” In Proceedings of the manila international sympo¬ 
sium on the development and management of energy resources, vol. 1, 
pp. 432-437. 1989. 

[4] Yan, M. Q., S. A. Sherif, G. T. Kridli, S. S. Lee, and M. M. Padki. 
’’Thermo-fluid analysis of solar chimneys.” In Industrial Applications of 
Fluid Mechanics-1991. Proceedings of the 112th ASME winter annual 
meeting, Atlanta, GA, pp. 125-130. 1991. 

[5] Von Backstrom, Theodor W., and Thomas P. Fluri. ’’Maximum fluid 
power condition in solar chimney power plantsan analytical approach.” 
Solar Energy 80, no. 11 (2006): 1417-1423. 


10 



[6] Schlaich, Jrg, Michael Robinson, and Frederick W. Schubert. The solar 
chimney: electricity from the sun. Geislingen, Germany: Axel Menges, 
(1995). 

[7] Putkaradze, Vakhtang, Peter Vorobieff, Andrea Mammoli, and Nima 
Fathi. ’’Inflatable free-standing flexible solar towers.” Solar Energy 98 
(2013): 85-98. 

[8] Fluri, T. P., and T. W. Von Backstrm. ’’Performance analysis of the 
power conversion unit of a solar chimney power plant.” Solar Energy 82, 
no. 11 (2008): 999-1008. 

[9] Peter Vorobieff, Andrea Mammoli, Nima Fathi, and Vakhtang 
Putkaradze. ’’Free-standing inflatable solar chimney: experiment and 
theory.” Bulletin of the American Physical Society 59 (2014). 

[10] Nima Fathi, Peter Vorobieff, Seyed Sobhan Aleyasin. ” V&V Exercise for 
a Solar Tower Power Plant.” ASME Verihcation and Validation Sympo¬ 
sium (2014). 

[11] Zhou, Xinping, Jiakuan Yang, Bo Xiao, Guoxiang Hon, and Fang Xing. 
’’Analysis of chimney height for solar chimney power plant.” Applied 
Thermal Engineering 29, no. 1 (2009): 178-185. 

[12] Koonsrisuk, Atit, and Tawit Ghitsomboon. ’’Theoretical turbine power 
yield in solar chimney power plants.” In Thermal Issues in Emerging 
Technologies Theory and Applications (ThETA), 2010 3rd International 
Gonference on, pp. 339-346. IEEE, 2010. 

[13] Xu, Guoliang, Tingzhen Ming, Yuan Pan, Fanlong Meng, and Gheng 
Zhou. ’’Numerical analysis on the performance of solar chimney power 
plant system.” Energy Gonversion and Management 52, no. 2 (2011): 
876-883. 

[14] Tayebi, Tahar, and Mahfoud Djezzar. ’’Numerical Analysis of Flows in 
a Solar Ghimney Power Plant with a Gurved Junction.” International 
Journal of Energy Science (2013). 

[15] Koonsrisuk, Atit, and Tawit Ghitsomboon. ’’Mathematical modeling of 
solar chimney power plants.” Energy 51 (2013): 314-322. 


11 



[16] Homayoon Daneshyar. ” One-dimensional compressible flow.” Oxford, 
Pergamon Press, Ltd., (1976). 

[17] Iain G Currie. ’’Fundamental mechanics of fluids”. McGraw-Hill, Inc., 
(1993). 


12 



